AN INVERSE PROBLEM FOR POINT INHOMOGENEITIES 
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Abstract. We study quantum scattering theory off n point inhomogeneities 
(n £ N) in three dimensions. The inhomogeneities (or generalized point inter- 
actions) positioned at {^i, . . . ,^n} C M.^ are modeled in terms of the (real) 
parameter family of self-adjoint extensions of —A| pjjo ^jj3^^^^ ^ jj inL^(R'^). 
The Green's function, the scattering solutions and the scattering amplitude for 
this model are explicitly computed in terms of elementary functions. Moreover, 
using the connection between fixed energy quantum scattering and acoustical 
scattering, the following inverse spectral result in acoustics is proved: The 
knowledge of the scattered field on a plane outside these point-like inhomo- 
geneities, with all inhomogeneities located on one side of the plane, uniquely 
determines the positions and boundary conditions associated with them. 



1. Introduction 

To describe the inverse problem solved in this paper in some detail we need a few 
preparations. Let Mj. = {x = a;2, 2:3) G M"^ | X3 ^ 0}, P = {x = {xi, X2, x^) G 
M'^ I ^3 = 0}, D C Mi^ a domain with smooth boundary and finitely many connected 
components, and v G L^{D), supp(t;) ^ D, v real- valued . Consider the fixed energy 
scattering problem 

{[- Vl -kl - klv{x)]u(x, y) = 5(x-y), x,y eM?,x^y, 
lim|j.|^oo y) — ifcou(a;, ?/)] — uniformly in directions uj = \x\~'^x 

and uniformly in y for y varying in compact sets, 

(1.1) 

with fco > (the wave number) a fixed positive constant. Here c{x) = [w(a;) + 1]^^/^ 
has the physical meaning of the wave velocity profile in the medium, v{x) is the 
inhomogeneity in the velocity profile, and u{x,y) represents the acoustic pressure 
generated by a point source at the point y G M.'^. 

The inverse problem (IP) associated with ( |l . l|) , more precisely, the inversion of 
the surface data u{x,y) for the velocity profile c{x), then can be formulated as 
follows: 



IP 1.1. Given the data {u{x,y)}x^yi£p at fixed ko > 0, determine v(x), x £ D. 

A solution of this inverse problem (i.e., uniqueness of v{x) and recovery of v{x) 
from the prescribed data) is described in Sects. III.6, IV. 2]. Numerical methods 
in connection with IP 1.1 are discussed in |ll|. Sect. V.3]. 
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Since u{x,y) in (1.1) can be identified with the Green's function at fixed energy 

fc2 > 0, 

G{klx,y) = {-A-klv-kl)-\x,y), x,yeR^x^y, (1.2) 

associated with the self-adjoint (Schrodinger-type) operator 

H = -A-klv, doin{H) = H^-^(W^) (1.3) 

in L^(]R'^), we can reformulate the inverse problem IP 1.1 in the following equivalent 
form: 

IP 1.1'. Given the data {G{kl,x,y)}x,y£p at a fixed energy k^ > 0, determine 

v{x), X <E D. 

For practical applications in connection with ultrasound mammography tests (as 
opposed to x-ray mammography) and in the area of material science in connection 
with the detection of cracks and cavities, it is of relevance to consider inhomo- 
geneities v{x) of the special form 

n 

v{x) = Y,v,{x - ^,), {6,...,Cn}cM3_, (1.4) 

i=i 

where Vj e L^{Dj), supp(z;j(- — C Dj, and Dj C are connected domains 
with smooth boundaries and sufficiently small diameters dj with respect to the wave 
length (i.e., (maxi<j<„ dj)fco 1). A numerical procedure recovering the S,j (and 
hence the approximate position of the small inhomogeneities) and the intensities 
of the inhomogeneities, defined by Vj — Jj^ dxvj{x — S.j), 1 < j < n, has recently 
been discussed in |]l^ . 

At this point we are in a position to describe the inverse problem considered 
in this paper. In view of the physical applications mentioned in connection with 



(1.4), we will now consider the idealized situation of inhomogeneities Vj{x — ^j) 
of point-like support at 1 < j < n. Intuitively, we want to solve the inverse 
problem 

'Vl -fco~fcoE"=i^^j(a;-Cj)]w(a;,y) = 5{x-y), a;,y €R^ x ^ y, 
lim|j.|^oo [g^w(x, y) — ikou{x, y)] — uniformly in directions cv — |a;|~"'^a; 

and uniformly in y for y varying in compact sets, 

(1.5) 



where formally 



- = aj'^(a^ - Cj), ^<j<n, (1.6) 



for some "coupling" constants a-,- € M, 1 < j < n. However, as is well-known, point- 



like inhomogeneities of the type (1.6) as potential coefficients in a Schrodinger-type 
operator in dimensions c? > 2 do not lead to an operator or quadratic form pertur- 
bation of the Laplacian -A, dom(-A) = H^''^{R''^) in L'^{R''-), where d denotes the 
corresponding space dimension. One possible way around this difficulty for d = 2 
and d = 3 is the introduction of an appropriate coupling constant renormalization 
procedure. This point of view is presented in detail in Ch. II. 1]. Alternatively to 
this renormalization procedure for d — 2,3, one can apply the theory of self-adjoint 
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extensions of closed symmetric densely defined linear operators in a Hilbert space 
to the operator 



-^lc,~(M3\H)' S = te,...,e„}cR-i (f.7) 

in L^{M!^). (Here the T denotes the operator closme of T and we refer to Remark 
for a brief discussion of the situation in different dimensions d G N.) In this paper 
we follow the latter approach and model the Laplacian —A perturbed by point- 
like perturbations of the type ~fco ~ ^j) self-adjoint extensions of 
— A|^o„^g3^„j, denoted by — A^.h, parametrized by the (real) parameter family 
of self-adjoint matrices 6* in C". 

Taking advantage of the equivalence of the inverse problems IP 1.1 and IP 1.1', we 
can now formulate the inverse problem associated with point-like inhomogeneities, 
as studied in this paper, in a precise manner as follows: 

IP 1.2. Prove that the data {Ge,s{kQ,x,y)}x,yep at fixed energy /cq > 0, uniquely 
determine S = {^i, . . . , C M'?_ and the self-adjoint n x n matrix 6 in C". 

Here Ge,s(kQ, x, y) denotes the Green's function associated with — A^.h, that is, 

Gff^siz, X, y) = i-Ae,s- zy\x,y), det(Pe,H(^)) ^ 0, x,y eR'^XE, x ^ y. 

(1.8) 

While IP 1.1 (resp., (IP 1.1') is concerned with uniqueness and reconstruction of 
v{x), X E D, IP 1.2, as studied in this paper, focuses on the unique determination 
of S and 9 by the data measured on the plane P. 

In Section ^ we present a detailed account of Krein's formula of self-adjoint 
extensions of closed symmetric operators in a Hilbert space, our principal tool in 
describing the ti^ (real) parameter family of self-adjoint extensions — A^.h of ( |1.7| ) in 
Section |^. In particular, we explicitly describe the Green's function, the scattering 
solutions, and the scattering amplitude associated with —Ag^s. in Section |^. The 
inverse problem IP 1.2 is solved in our final Section ^j. 



2. Krein's Formula for Self-adjoint Extensions 

In this section we recall Krein's formula, which describes the resolvent difference 
of two self-adjoint extensions Ai and A2 of a densely defined closed symmetric 
linear operator A with deficiency indices {n,n), n G N. (Reference treats this 
topic in the general case where n S NU {00}. Here we restrict ourselves to the case 
71 < 00.) We start with the basic setup following 

Let H he a separable complex Hilbert space, A : dom(^) — > H, dom(^) = 
Ti a densely defined closed symmetric linear operator in Ti. with finite and equal 
deficiency indices def(A) = (r, r), r e N. Let Ai, i = 1,2, be two distinct self- 
adjoint extensions of A and denote by A the maximal common part of Ai and A2, 
that is, A is the largest closed extension of A with dom(^) — dom(^i) ndom(^2)- 
In this case one calls Ai and A2 relatively prime with respect to A. Let < p < r — 1 
be the maximal number of elements in dom(^) = dom(^i) n dom(yl2) which are 
linearly independent modulo dom(A). Then A has deficiency indices def(^) = 
{n, n), n — r ~ p. Next, denote by ker(A* — z), z E C\M the deficiency subspaces 
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of A and define 

VFi,,,,„ =I+iz- zo)iAi - = {Ai - zo){Ai - z)-\ z, zo G p(Ai), (2.f ) 

wfiere / denotes the identity operator in Ti. and p(T) abbreviates the resolvent set 
of T. One verifies 

^■1,^0,^,^1,^1,^2 = W^i,zo,z2 7 zo,zi,Z2 e p(Ai) (2.2) 

and 

Wi,.,.o ker(A* - zo) = ker(A* - z). (2.3) 

Let {uj{i)}i<j<n be an orthonormal basis for ker(A* — i) and define 

Ujiz) = Vl^i,^,»Uj(i) = {Ai - i){Ai - zy^Ujii), l<j<n, p{Ai). (2.4) 

Then {uj{z)}i<j<n is a basis for ker(A* — z), z G p(^i) and since = 
(Ai — i)(v4i + i)^^ is the unitary Cayley transform of Ai, {uj{—i)}i<j<n is in fact 
an orthonormal basis for ker(A* + i). 

The basic result on Krein's formula, as presented by Akhiezer and Glazman 
Sect. 84, then reads as follows. 

Theorem 2.1. (Krein's formula, [§ Sect. 84].) 

There exists a Pi.2{z) — (Pi.2{z)j,j') ,^ ■ G Af„(C), z G p{A2) D p{Ai), such 
that 

det(Pi,2(z)) 7^ 0, z G p{A2) n p{Ai), (2.5) 

PiA^r' - A,2(2o)"' zo)((w,(z),u,v(zo)))^^^.^.,^^ , z,zo G p(Ai), (2.6) 

Im(Fl,2(^)-^) = -/„, (2.7) 

n 

(A2-z)-i = (Ai-z)-i+ ^ Pi^2{z)j,j'{u,,{z),-)u,{z), zep{A2)np{Ai). 

(2.8) 

Here Im(r) = {T-T*)/{2i) and Re(r) = {T + T*)/2 denote the imaginary and 
real parts of the matrix T, respectively. 

We note that ^1.2(2:)^^ extends by continuity from z G p{A2) fl p{Ai) to all of 
piAi) since the ri ght- hand side of ( ^.6D is continuous for z G piAi). The normal- 
ization condition (|2.7| ) is not mentioned in |Q but it trivially follows from (^^) and 
the fact 

=^«', l<i,/<n (2.9) 
(where 5jji denotes Kronecker's symbol) and from 

Pi.2{z) = PiAz), z G p{Ai) n p(A2). (2.10) 



Taking z = zq in ( ^.6[ ) shows that —Pi, 2 (2) ^ and hence ^'1.2(2:) is a matrix- valued 
Herglotz function, that is, 

Im(Pi,2(z)) > 0, z G C+. (2.11) 



Strict positive definiteness in (2.11) follows from the fact that {uj(z)}i<j<„ are 
linearly independent for z G C+ and hence {{uj{z),Uj'{z))i<cjji<n > 0. 

Next we turn to the connection between Pi^2{z) and von Neumann's parametriza- 
tion of self-adjoint extensions of A as discussed in detail in Q. Due to (2.6), 
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-Pi, 2(2) ^ is determined for all z G p{Ax) in terms of Pi^iii) ^, (^1 — z) ^ and 
{Mj(i)}i<j<„, 

Pi,2(^)-' = Pi,2(^)-' -{z- i)/„ - (1 + (Ai - zY^u,,{i)))^^^^^,^^^, 

zep(Ai). (2.12) 

Hence it suffices to focus on 

Fi,2«-' - Re {Pi.i(i)-^) - *^n- (2-13) 

Let 

:ker(A* -i) ->ker(A*+i), ^=1,2, (2.14) 

be the linear isometric isomorphisms that parameterize Ai according to von Neu- 
mann's formula 

A,(/ + (/ + W,)ii+) =^/ + ^(/-^^,)u+, 

dom(A/;) = (/ + W£)u+) e dom(A*) I5 e dom(A), u+ G ker(yl* - i)}, 

£=1,2. (2.15) 

Next, denote by Ui = (C^£j,i')i<ij'<n G Af„(C), £ = 1,2 the unitary matrix rep- 
resentation oilAi with respect to the bases {Mj(«)}i<j<ri and {ui j(— i)}i<j<„ of 
ker(A* — i) and ker(yl* -I- i) respectively, that is, 

n 

UiUj{i) = ^ Ui^y.jU^.y^-i), l<j<n, £ = 1,2. (2.16) 

Lemma 2.2. (0.) 

(i) t/i = 

(ii) -1 ^ spec(C/2). 

(iii) Ui, £=1,2 and Pi,2{i) are connected by 

Pl,2(*) = \{In + = ^(t/2"' - C/r'). (2.17) 

Here spec(r) denotes the spectrum of T. 
Next, writing 

[72 = exp(i02), ^2=^2 (2.18) 

for the matrix representation of U2 with respect to the bases {uj(i)}i<j<„ and 
{uij{—i)}i<j<n of ker(A* — i) and ker(A* -|- i), one verifies 

Re(Pi,2(i)"') = tan(02/2). (2.19) 

Introducing the matrix- valued Herglotz function Mi{z) associated with Ai (cf. Q, 
0)by 

Mi(z) = z/„ + (1 + z')((^iJ(^), (Ai - z)"'"/ («))i<j j'<n) , 2 e p{Ai), (2.20) 

Pi.2(-z) in Krein's formula (2. J) then can be rewritten as 

Pl,2(z) = (tan(^^2/2)-A/l(z))-l 

= ( tan(6'2/2) - z/„ - (1 + z2)((Mj(i), (Ai - z)"^Wj'(i))i<jj'<„))"\ 

0ep(Ai). (2.21) 
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We emphasize that 



{((1/2) +TO)7r}„ez i spec(6'2) 



(2.22) 



according to Lemma 2.2 (ii), due to our hypothesis that A\ and A2 are relatively 
prime with respect to A. 

For subsequent purposes it is useful to introduce the self-adjoint operator ■&2 G 
i3(ker(A* — i)) defined through its matrix representation Q2 with respect to the basis 

{wj(i)}i<j<„, that is, 



(2.23) 



The discussion of Krein's formula thus far dealt exclusively with the orthonormal 
bases {uj{i)}i<j<n and {ui j (— «)}i<j<„ of ker(v4* — i) and ker(A* + i) following 
our discussion in Q and Appendix B]. In the remainder of this paper, however, 
it will be be more convenient to discuss matrix represe nta tions of Mi{z) and 1^2 
with respect to a natural (cf. the comment following (3^)), but not necessarily 
orthogonal basis. Hence we b riefl y discuss the effect of a change of basis in con- 
nection with Krein's formula (2^). Let {ujii)}i<j<n be another (not necessarily 
orthogonal basis) of kcr{A* — i) and define 



Uj{z) = {Ai - i){Ai - z) ^Uj{i), 1 < j < n, z G p{Ai), 

n 

UiUj{i) = ^ Uij'jUj>{-i), 1 <j <n, ^=1,2, 

C/2 = exp(i02), e*2^92. 
In addition, one verifies 

Ui - -/„ 



(2.24) 
(2.25) 

(2.26) 
(2.27) 



as in Lemma 2.2 (i). Krein's formula (2.8) then can be rewritten in the form 



iA2 



'+ PiAzhri^r{z),-)ujiz), z e p{A2) n p{Ai), 



(2.28) 

where 

Pi^2(z) = (tan(^2/2)-Mi(z))-i 

= ( tan(^2/2) - zin - (1 + Z^){{Uj{i), {Ai - 2;)"^Wj'(i))l<ij'<n))~\ 

zep(Ai), (2.29) 

and (cf. d^H)) 

= ^2Uf{t)), 1 < j,f < n. (2.30) 

The proof of (|2.28|) -( ^.3C| ) is based on the following elementary result. 

Lemma 2.3. Let Tin, n G N be an n- dimensional complex Hilbert space, T 6 B{7i) 
a bounded linear operator in TL with G B{7i). Assume that {ipj}i<j<n and 
{'>pj}i<j<n o,re {not necessarily orthogonal) bases in Ti. Then 

n 

^ ((('0^r'0m)l<£,m<„) ) (V'i' , OV'i 
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= X! { {i^i^T'^rn)l<£^m<n) ^) ^. ("0/ , • ■ (2-31) 

3. The Direct Scattering Problem for Generalized 
Point Interactions 

In the principal part of this section we apply the abstract framework surrounding 
Krein's formula ( 2.2^ ) to the concrete situation of n generalized point interactions 
in M.^ . At the end we derive the corresponding quantum mechanical scattering 
formalism, including explicit expressions for the scattering wave functions and the 
scattering amplitude. 

In order to apply the results of Section we now make a series of identifications: 

H = L'^{m.^), (3.1) 



Cs^{R^\{iu-.iJ ' ^^1' ■ • • '^"} ^ for ] + j\ (3.2) 

ker(A* -i) = span{{tj (i, x) = Go{i,x - Cj)}i<j<ri, (3.3) 



N-i/ exp{iz^^^\x ~ y\) 



Go(z,x-y) = (-A-z)-^(a;,y) = ^:^^^^, (3.4) 

4:Tt\x — y\ 

z e C\M, Im(zi/2-) > 0, e R^ a; ^ y, 

Ai = -A, dom(-A) = B^^'^i^), (3.5) 

u,{z) = (-A - z)(-A - z)-^u,(i) = Go{z, ■ - 0), (3.6) 

z e C\R, Im(zi/2) > 0, \ <3<n. 

In particular, a comparison of ( |3.3D and 

ker(A* — z) = span{uj(z, x) = Go(2, a; — ^j)}i<j<„, z £ C\R (3.7) 

shows that \uj{z)\\<ij<_n is a n atural (though, not orthogonal) basis of ker(^* — z). 
We note that the fact (^^) can be found, for instance, in jH Sect. II. 1.1] and 

Straightforward computations using 



(±i)l/2 = 2-l/2(±l+j)^ i(i)l/2 =j(_j)l/2 (3,8) 

and the first resolvent equation 

(-A - zi)-n-A - z^)-^ = (zi - Z2)-M(-A - zi)-i - (-A - z2)-i], (3.9) 

zi, Z2 e p(-A) 

repeatedly, then yield the following results. 

Lemma 3.1. Let z G p(— A) and j, j' e {1, . . . , n}. Then 

(uj(0,u/(«)) = < , ^r- r f c ^^ • v ^^'^^^ 

(u,(z), (-A - z)-^u,{i)) = (47r)-i(z2 + l)-i[izi/2 - i^-if^^ - 2~^'^{z + i)], 

(3.11) 

({t,(z), (-A - z)-%,{i)) - (z2 + l)-MGo(z,^, - 0') - Go(-*,e, - ^r)] 

-(z-z)-ilm(Go(z,0-ey)), jV/, (3.12) 
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where 

Im(Go(z,C, - ^r)) = [Go(^,^ - ^r) - Go{z,^j - ^r)]/i2^). (3.14) 

Given these preliminaries, one can now describe the (real) parameter family 
of all self-adjoint extensions of A, rela tiv ely p rime to Ai — —A with respect to A, 
by appealing to Krein's form ula ( ^.28| ), ( ^.29| ) as follows. (In passing we note that 



Ai — —A, as defined in (3^), is the Friedrichs extension of A.) One defines 

S = Ui,...,^„}cM3 (3.15) 

and denotes by 



a self-adjoint n x n matrix in C". Combining Krein's formula ( 2.2S| ), ( ^.29| ) with 



(3.11) and (3.12) then yields the principal result of this section. 



Theorem 3.2. Let z G C\IR.. Then the ri^ [real) parameter family of all self- 
adjoint extensions — Ag.s of — Al^^^j^g^..^ , relatively prime to —A with respect to 

— A|^^^jj3^„^ , can be parametrized by all self- adjoint n x n matrices 9 in C" with 

{{{l/2)+m)n},nezispec{9). (3.17) 
An explicit representation for — A^^^ is provided by 

n 

(-A,,H - z)-^ = (-A - z)-i + Pe,E{zhAGo{z, ■ - 0'), ■)Go{z, ■ - C,), 

(3.18) 



whe 



PeMz)-')^^^, (3.19) 



and 



-(47r)-izzi/2 _ (4^)-i2-i/2 _^ (tan(6i/2)jj, 
-Go(^,0 - ^r) + Re(Go(*,0 - ^f)) + (tan(0/2)),-y, 3 + 3', 

Re(Go(^,0 - 0')) = [Go{z,^, - C/) + Go(z,^, - Cj')]/2. (3.20) 



Remark 3.3. (i) Whenever 9 in (3.18|) has an eigenvalue ((1/2) + mo)7r for some 



rriQ £ Z, Pe^z) becomes a singular matrix, det(Pe^H(z)) 0, z G C\M. In this case 
at least one point is removed from S and one effectively considers self-adjoint 
extensions of A = — A| „^ ,„3, ,^ , , , , parametrized in terms of (n — 1) x (n — 1) 
(or less) dimensional self-adjoint matrices 9. In particular, the Friedrichs extension 



of A = —A goo^]j^ 3^~-j, given by Ai = —A, formally corresponds to the extreme case 



9 = 7r/„ in (|1|), (|3l9| ) 



(ii) It seems appropriate to call the n^-parameter family — Ag^s defined by (3.18), 
( 3.19| ) the generalized point interaction Hamiltonian, distinguishing it from the usu- 



ally considered n-parameter family of (local) point interactions. In fact, introducing 
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a — (ai, . . . , an) G M", the standard n-parameter family of self-adjoint extensions 
—Aa,s of — A|p^^jj3^„^ emerges as a special case of ( 3.18 ), ( 3.19 ) by choosing 

tan(^/2) = {{a, + {4n)-^2~^/')S,^,, - Re(Go(i,C, - C/)))i<,,,,<„, (3.21) 
where 

Go(z,x) = |^''^''"^ , (3.22) 

^ ^ jo if a; = (or if n = 1), ^ ' 



andRe(Go(^,Cj -Cj')) = [Go(^, C - 6') + Go(z, - Ci')]/2- Insertion of (|^ into 
( p^ ), ( |37l9| ) then yields 

(-A„,H - z)-^ - (-A - z)-i + ^ (r^^H(^)-^)^. ^.,(Go(2, • - ^,0, ■)Go{z, ■ - ^,), 

(3.23) 

where 

ra,H(^) - {{-{inrhz'^^ + a,)S,,,> - G^{z,i, - )) i<,- ,,<„, (3.24) 

in accordance with ^, p. 113]. While most efforts in connection with finitely many 
point interactions focus on the n-parameter family — Aq^h (cf. the detailed discus- 
sion in H, Ch. II. 1] and the references therein), the general ri^-parameter family of 
generalized point interactions has been discussed by Dabrowski and Grosse js) in 
1985. The treatment in ||] also combines Krein's resolvent formula with von Neu- 
mann's parametrization of self-adjoint extensions, but is somewhat less detailed 
than our present approach. (In particular, their matrix 5'(z, zq), and hence their 
M(z), are not explicitly computed in section II of ||], although these quantities can 
be inferred from the scaling limit approach in section IV via their formula (4.18).) 

Finally, we briefly discuss stationary quantum scattering theory foll owin g the 
lines of Sect. II. 1.5] and ||l^. Given the resolvent kernel of — Ag 5; in ( |3.18 ), one 
computes 

lim lim 47r|?/|e-*('=+*")l«l(-Ae,H - (fc + ie)2)"^(a;,y) 

e^O \y\ — >oo 

n 

= e''"^-'=+ Pe,E{k^)j,j'e'''"-^^'Go{k^,x-^/) (3.25) 

= ^e^x, k, uj), keR, det(Pe,H(fc^)) ^0, iv e S^, x e R^\E . 
Moreover, since 

{-A^){x,k,Lu) ^ k^-^{x,k,uj), xeR^\E (3.26) 

in the distributional sense as well as pointwise, '^g,-s{kuj, x), k €M., det(Pe^H(fc^)) 7^ 
0, UJ € S"^, X G ]R"^\S, represent the generalized eigenfunctions, that is, the quantum 
scattering wave functions associated with — A^ h. 

The corresponding quantum scattering amplitude ^6.3(0;', tj, fc) is then com- 
puted as follows, 

Ae,H(w',w,/c) = lim \x\e-'''\'=\[i'e^s{x,k,uj)-e''"^-''] 

\x\ — 'OO 
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n 

= (47r)-i J2 Pe,B.{k^)j,r(^'^^""^'''^''^'\ (3-27) 

A: e M, det(Pe,H(fc^)) 7^ 0, t^,t^' £ S^. 
The corresponding scattering matrix Sg,-E.(k) in L'^{S'^) is then given by 

Sem-I+^, E ^«,H(fc^W(e-'^-«^"(-\>-^'=«-(-), (3.28) 
jj''=i 

fc e R, det(Pe,H(fc^)) 7^ 0. 

Remark 3.4. Since Sg,-s{k) is unitary in L'^{S'^) (this either follows from abstract 
methods since —A and — Ae^^ are self-adjoint and the second term on the right- 
hand side of ( 3.18| ) is of rank n and hence a trace class operator, or directly from 
( ^.19| ) and ( 3.2§| )), the scattering amplitude Ae;E{i-^' ^^^Tk) (the integral kernel of 
Se^sik) — /) automatically satisfies the (generalized) optical theorem, that is, 



Im(Ae,H(w',w,fc)) = (47r)-ifc / doj" Ag^sii^" ,uj,k)Ae,s{u;",uj',k) , (3.29) 

/c e R, det(P0,H(A:^)) ^ 0, uj,uj' e 5'^. 
On the other hand, reciprocity of the scattering amplitude Ag 5(0^', uj, fc), defined 

by 

^e,H(w', w, fc) = Ag^si-i^, -uj', k), fee E, det(Pe,H(fc^)) 7^ 0, w, to' e S'^ (3.30) 
is satisfied if and only if 

6** = 0, (3.31) 

where T* denotes the transpose of the matrix T. Together with the requirement of 
self-adjointness of 9, 0* = 0, this yields an n{n + l)/2 (real) parameter family of 
operators —Ag^s satisfying 0^9* = 0* . (The number of real elements above and 
on the diagonal of equals X]j=i J — "■('^ + l)/2.) 

Similarly, the reality constraint on Ag^s.i'^'i^, k), that is, the requirement 



Ag^sito', UJ, k) = Ag,;E{uj', , ~k), k E R, det{Pg^s{k^)) ^ 0, UJ, Lj' e 52 (3.32) 

is satisfied if and only if is a real matrix, 

hr='^^ l<J,/<". (3.33) 

Together with self-adjointness of this again results in — 0* — 0* and hence is 
equivalent to the reciprocity requirement. (For background material on properties 
of the scattering amplitude, such as the optical theorem, reciprocity, and reality, 
we refer to Sect. 1.4] for obstacle scattering and [|l3[ Sect. 3.6] in the context of 
potential scattering.) 

It is interesting to observe that these natural requirements on the scattering 
amplitude, such as the optical theorem, reciprocity, and reality, are satisfied for an 
n{n + l)/2-parameter family of generalized point interactions (though, not for the 
full n^-parameter family) and hence for a larger family than the usually considered 
n-parameter family of (local) point interactions —Aa^s.- 

We conclude this section with the following remark on space dimensions other 
than three (the interested reader can find many more details in ]S]). 
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Remark 3.5. All results of this section immediately extend to the case of two space 
dimensions replacing the Green's function (3.4) of the three-dimensional Laplacian 
by the corresponding two-dimensional Green's function 

Go(z, x-y)^i-A- z)-\x, y) = {i/A)hI^\z^/^\x - y\), (3.34) 

Im(z^''^) > 0, a;, y e M^, a; / y. 

Here H^pi^-) denotes the Hankel function of order zer o an d first kind (c f. O 
Sect. 9.1]). There are only m inor changes required in ( 3.1C ), ( 3.11 ), and ( |3.12|) 
due to the ln(z)-behavior of ( 3.34 ) as z ^ 0. Analogous results apply to the 
one-dimensional case using 

Go{z,x- y) = i-A ~ z)-\x, y) - {i/2)z-''\'^"'"\--y\ (3.35) 

Im(z^/^) > 0, e M. 

The one-dimensional case, however, differs from the two and three-dimensional cases 

since B — [(^oo^jj^j-jj^ ^ j-, now has deficiency indices (2rt, 2n) (as opposed 

to {n,n) in two and three dimensions, cf. (|3.3|) ). Consequently, B admits a 4n^ 
(real) parameter family of self-adjoint extensions and hence additional types of 
(generalized) point interactions in dimension one. The proper definition of A with 

ker(A* - z) = span{G'o(2, • - Cj)}i<j<n (3.36) 

in one dimension is given by 

dom{A) = {geH^-\R)\g{Cj)=0,l<j<n}. (3.37) 

Finally, since 



-^lc„-(E^\{?i,.. ..«„}) ^ -^l//2,2(R.) for d > 4 (3.38) 

(i.e., — A|p^^jj^^^^^ ^ is essentially self-adjoint for d > 4), there are no (gener- 
alized) point interactions in four dimensions or higher. 

4. A Uniqueness Result 

Given the preparations in Section ^, the principal purpose of our final Section |^ 
is to provide a solution of the inverse problem IP 1.2 formulated in Section |l|. More 
precisely, we will prove the following uniqueness result (we freely use the notation 
established in Sections |l|-^ throughout this section). 

Theorem 4.1. Let fco > and assume that dct{Pg^s{ko)) ^ 0- Then the data 
{Gg^-E.{kQ, x,y)}x,y(^p uniquely determine S = {Cij---i^n} ^'^'^ ^^V' 

adjoint n x n matrix 9 in C" . 

Proof. Given the data Gg^E{kQ,x,y) for all x,y e P, x ^ y, det(Pe^H(fco)) 7^ 0, the 
Poisson-type formula, 

f d 

w{s, y) = J Ge^s{kl, (cr, 0), ?/)— [Go(fco, (ct, erg) - s) 



-Go(fco',(^,-^3)-s)]|^^^„, (4.1) 
s e K^, yeP 
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(here a — (cri, (72) G M ) yields the solution of the problem 



(- -fc§)w(a;,y) = 0, xe 



lim|2.|^oo y) — ikQw{x,y)'\ — uniformly in directions lu — \x\^^x 

and uniformly in y for y varying in compact sets 



(4.2) 



Ge,E(fco. 2;, y) for aU a; e U P, y G P, a; 7^ y. (4.3) 



for each fixed y E P. In particular, w{x, y) in (4.1) represents 

By symmetry of the Green's function Ge,H(fco, a;, y) with respect to a; and y, 

GeMko,x,y) =Ge,Eikl,y,x), x,y eR^, x^y, (4.4) 

we also determined 

Ge,s{kl,x, y) for all a; e P, y e U P, a; ^ y. (4.5) 

Moreover, using Ge^~{kl, (ct, 0), y) with y g U P (instead of y G P) in ( |4l| ) then 
determines 

Ge^H (fco ' 2;, y ) for all X, y G R^ U P, X ^ y . (4.6) 

In other words, we managed to lift the data from P to R'^ U P. 
Next, the explicit formula ( |3.18| ) for (— Ag s — z)~^ yields 

n 

Ge,~.{z,x,y) = Go{z,x,y) + ^ Pe^3{z)j,j'Go{z,x - (j)G){z,y - (4.7) 

det(Pe,H(z)) ^ 0, a;, y G R^S, a; / y, 



with P6(,h(z) defined in (3.19). Hence one concludes 

(- v' -2)Ge,H(^, X, y) - 0, x, y G R^S, x ^ y. (4.8) 
Thus, the data {Ge^siko, x,y)}x,y£P, det{Pg ^^{kQ)) ^ uniquely determine 

G9,s{klx,y) for aU x,y G R^S, x^y (4.9) 

by the unique continuation property ||, Sect. 17.2] applied to (4^). The singularity 
structure of (4.7) then determines Ci: ■ ■ • i ^-nd hence S. Similarly, taking x — > 
and y — > ^j' independently, determines Pd,s{kQ)jj>, 1 < < n, and hence 0. 
Thus, Theorem 4.1 is proved. □ 
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